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(54) METHODE POUR SIMPLIFIER LA MODELISATION D'UN MILIEU GEOLOGIQUE POREUX TRAVERSE PAR UN 
w RESEAU IRREGULIER DE FRACTURES. 



(57) Methode permet une modelisation simplifiee d'un mi- 
lieu geologique heterogene poreux (tel qu'un reservoir tra- 
verse par un reseau irregulier de fissures par exempie) 
sous la forme d'un milieu transpose ou equivalent de ma- 
niere a ce que te milieu transpose soit equivalent au milieu 
d'origine. relativement a un type determine de fonction de 
transfer! physique (connu pour le milieu transpose). Elle 
comporte: 

<a) la formation d'une image en au moins deux dimen- 
sions du milieu geologique d'origine sous la forme d'une 
l serie de pixels et I'association a chaque pixel de la serie 
d'une valeur initiale particuliere pour cette fonction. 
b) la determination par etapes de valeurs a attribuer pour 
LO la fonction de transfert physique a chaque pixel de la serie 
rt) (telle que la distance minimale separant le pixel de la fis- 
sure la plus proche) par reference a des valeurs de la fonc- 
fs» tion attributes a des pixels voisins de I'image, et 
jrt c) la determination d'une propriete physique du milieu 
■s^ transpose ou equivalent en identifiant ies valeurs de la 
1 fonction de transfert connue pour !e milieu transpose avec 
£^ les valeurs de la fonction de transfert du milieu d'origine 
determinees par etapes. La fonction de transfert physique 
fy* peut representer des variations entre dtfferentes parties du 
milieu geologique. par exempie des distances ou des 
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transmissivites ou des transferts thermiques (entre un re- 
servoir et un puits traversant te reservoir par exempie). etc. 

Application a la determination d'un milieu transpose of- 
frant !a meme recuperation d'un fluide au cours d'un pro- 
cessus d'imbibition capillaire que le milieu reel. 
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^invention concerns une meihode desunee a simpliiier !a moderation dun milieu 
geoiogique poreux traverse par un reseau irregulier de fissures qui simplitle la m.se en 
relation de modeles de caracterisation de reservoirs fissures et de modeles double 
porosite. La methode peut etre mise en oeuvre par exemple dans le domaine de la 
production petroliere par des ingenieurs de gisement en vue d'obtenir des predictions 
decoulement fiables. 

Les reservoirs fissures constituent un type extreme de reservoirs heterogenes 
comportant deux milieux contrasts, un milieu matriciel comenam la plus grande part de 
I'huile en place et presentant une faible perrneabilite. et un milieu fissure representant 
moins de I % de I'huile en place et hauternent conducteur. Le milieu fissure lui-meme 
peut etre complexe. avec differems ensembles de fissures characterises par leur densite, 
longueur, orientation, inclinaison et ouverture respectives. Les images en 3D de 
reservoirs fissures ne peuvent pas etre utilisees directement sous forme de donnees 
d'entree de simulation de reservoir. La representation dun reseau de fissures dans des 
simuiateurs decoulement a ete iongtemps considered comme irrealiste car la 
configuration du reseau est en partie inconnue et a cause des limitations numeriques liees 
a la juxtaposition de nombreuses cellules presentant des dimensions et des proprietes 
extremement contrasted. C'est pourquoi la modelisation simplifiee mais realiste de tels 
milieux presente un grand interet pour les ingenieurs de gisement. 

L" approche double porosite" telle quelle est enseignee, par exemple. par Warren J.E. 
et al dans 'The Behavior of Naturally Fractured Reservoirs". SPE Journal (septembre 
1963 ). 245-255. est connue dans Tart pour interpreter le comportement d un ecoulemem 
monophasique observe en testant un reservoir fissure. Selon ce modele de base, tout 
volume elemental du reservoir fissure est modelise sous la forme dun ensemble de 
blocs parallelepipediques identities Nmite par un systeme orthogonal de fissures 
uniformcs continues oricntees dans la direction de Tun des iro.s principaux sens 
decoulement. Lecoulcment des fluides a lechelle du reservoir s'efTeui.e a (ravers le 
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milieu de fissure seulement et des echanges de fluidcs surviennent tocalemem entre les 
fissures el Ies blocs malriciels. 

De nomhrcux simutateuis de reservoirs fissures utii ete developpes en uiilisum un lei 
modele. a\ec des ameliorations ^peci fit] Lies relatives a la mode iisaiion des echanges de 
5 fluides entre mairiee et fissure regis par des forces capillaires. gravitationnellcs. 
\ isqueu.se s et des mecanismes eompositionnels. ainsi que des echanges matrice-niatrice 
(simulateurs double permeabilite double porositc). De nombreux exemples de 
techniques de Part anterieur sont cites dans Ies references indiquees ci-apres. 

- Thomas. L.K. et al : "Fractured Reservoir Simulation". SPE Journal (fevrier 1983). 
10 42-54; 

- Quandalle. P. et al : "Typical Features of a New Multipurpose Reservoir 
Simulator", SPE 16007 presente au 9eme Symposium sur la simulation de reservoirs de 
San Antonio. Texas. 1-4 fevrier 1987; 

- Coats. K.H. : "Implicit Compositional Simulation of Single-Porosity and Dual- 
ly Porosity Reservoirs". SPE 18427 presente au Symposium SPE sur la simulation de 

reservoirs de Houston. Texas. 6-8 fevrier 1989. 

L'un des problcmes que rencontrent les ingenieurs de gisement consiste a doter ce 
modele de base de parametres en vue d'obtenir des predictions d ecouiemem fiables. En 
particulier. les permeabilites de fissures equivaientes ct la Lai 1 le des blocs matricicls 

20 doivent etre connues pour chaque cellule du simulateur decoulement. Mors que la 
permeabilite matnciellc pent etre evaluee a partir de carottes. les permeabilites du reseau 
de fissures contenu dans la cellule, e'est-a-dire les permeabilites de fissures 
equivaientes, ne peuvent pas etre evaluees simplcment et necessitent !a prise en compte 
de la geometric et des proprietes du reseau de fissures reel. Une methodc permettant de 

25 determiner les permeabilites de fissures equivaientes d'un reseau de fissures est 
presentee dans la demande de brevet parallele EN. 96/16330. 

11 ex isle une procedure de reference connuc permettant de determiner les dimensions 
a.h de chaque blue d'une section travcrsee par un mailkige regulier de fissures i : eq qui 
est cquivalcnte a la section d'un milieu multicouchc nature! fissure tra\cr^e par un 
*' ' reseau de fissures FN ie long d'un plan de reference parallele aux couches (ec plan elant 



2757957 



habiiuellemem horizontal ou sensiblement horizontal). Puur chaque couche du volume de 
roche fissuree etudie {figure I ). les dimensions "horizontales" a, h des blocs do la section 
eqimalente sont determiners de nianiere iterative au mo\en de calculs et par 
comparaison des fonciions de recuperation d'huile en Ibnciion du temps Rm et Reqm 
respective men t dans !u section reelle RE du volume de roc he fissuree etudie et dans !a 
section EQ decoupee en "morccuux de sucre" (sugar lumps) de taille ideniique equivalent 
a la distribution des blocs reels. Cette methode convent ionnelle necessite un simulaieur 
decoulement polyphasique simple porosite discretisant les blocs matriciels et les fissures 
de nianiere a pouvoir comparer les courbes de recuperation. Une telle procedure est 
extremement couteu.se dans la mesure ou la discretisation de la section reelle pent 
comporter un nombre tres eleve de cellules. En fait, la forme reelle des blocs doit etre 
representee au moyen de cellules de fissures fines le long des limites de chaque bloc. La 
mat rice doit egalement etre discretisee avec un nombre suffisant de cellules en vue 
d'obtenir une fonction de transfer! d'imbibition bloc-fissure precise. 

On peut citer differentes techniques de 1'art anterieur mises en oeuvre dans ce 
domaine. parexemple : 

- Bourbiaux. B. et al : "Experimental Study of Cocurrent and Countercurrent Flows in 
Natural Porous Media". SPE Reservoir Engineering (aout 1990). 361-368. 

- Cuiec. L. et al : "Oil Recovery by Imbibition in Low-Permeability Chalk", SPE 
Formation Evaluation fseptembre 1994). 200-208. 

Cependant. les caracteristiqu.es specifiques d'imbibition n'ont encore jamais ete 
utilisees pour trouver les dimensions du bloc equivalent dans les modeles double 
porosite. Ainsi. les ingenieurs de gisement manquent dun outil systemattque permettant 
de calculer les dimensions des blocs parallelepipediques equivalents a des ecoulements 
polyphasiques pour determiner la distribution reelle des blocs dans chaque zone de 
reservoir fissuree. 

On connait egalement dans Tart anterieur des techniques permettant d'intesrer des 
donnees de llssuraiion naturelle dans des modeles de reservoirs fissures. Les donnees de 
rissuration sont notamment des donnees de nature gcometnque et elles component des 
niesures de densite. dc longueur, d'azimui et d'inclinaison des plans de fissure soil 
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observes sue des affleurements. duns des galeries dc mines, sur des carones. soil deduits 
a part if de diagraphies. Different?* ensembles de fissures peuvent el re disi indues et 
caracterises par differenies distributions statist iques des caraeteres des fissures. Les 
configurations des fissures une fois caracterisees. on peut realiser ties reseau x 
ruinieriqu.es de ces ensembles dc fissures en utilisant un proeede s;ochast iqiie re^pectant 
ies distributions Matistiques de> parametres de fissures. De eels procedes sont presented 
par exemple dans les brevets FR-A-2.725.S 1 4: 2.725.794 ou 2.733.07? du demandeur. 

La- methode selon 1'invention permet une moderation simplified dun milieu 
geologique d'origine heterogene poreux. (tel qu'un reservoir traverse par un reseau 
irreaulier de fissures par exemple j sous la forme d'un milieu transpose ou equivalent- dc 
telle maniere que le milieu transpose soil equivalent au milieu d'origine relativement a un 
type determine de fonction de transfert physique connue pour le milieu transpose. La 
methode comporte : 

- la formation dune image a au moins deux dimensions du milieu geologique sous la 
forme dune serie de pixels et Vassociation a chaque pixel de la serie d'une valeur initiale 
particuliere pour ladite fonction, 

- la determination par etapes de la vaieur a attribuer pour la fonction de transfert 
physique a chaque pixel de ladite serie. en reference a des valeurs de la fonction 
attributes a des pixels voisins de 1'image. et 

- la determination d'une propriete physique du milieu transpose ou equivalent par 
identification des valeurs de la fonction de transfert connue pour le milieu transpose 
(simplifie) avec la valeur de la fonction de transfer! determined par etapes pour le milieu 
d'origine. 

La fonction de transfert physique peut representer des variations entre differenies 
parties du milieu geologique. par exemple des variations de distances, de transmissivit^s 
ou de chaleur dels que des transfers thermiques entre un reservoir et un puits truversant 
ce reservoir}, etc. 

La methode peut ctrc appliquce par exemple pour determiner, a pariir d'une image 
dun milieu geologique poreux traverse par un reseau irrcgulier dc llssurcs. un milieu 
transpose comport ant un ensemble de blocs disposes de maniere regulierc et separes par 
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un maillage regulier de figures, lequel milieu transpose donnc sensibiement une mo me 
recuperation de fluide au cours dun processus d'imbibition eapillaire que le milieu reel. 
Dans ce cas. la methode com port e : 

- la tormation dune image a au moins deux dimensions du milieu red sous la forme 
5 d'une serie de pixels. 

- la determination, pour chaque pixel, de la distance minimale separant le pixel de la 
figure la plus proche. 

- la formation d une distribution de nombres de pixels en fonction des distances 
minimales aux fissures et la determination, a partir de cette distribution, de la fonction de 

10 recuperation (R) dudit ensemble de blocs, et 

- la determination de dimensions (a, b) des blocs reguliers equivalents de 1'ensemble a 
partir de la fonction de recuperation (R) et de la fonction de recuperation (Req) de 
(equivalent ten utilisant par exernpie une procedure d'identification desdites fonctions de 
recuperation). 

15 Avec la methode definie ci-dessus utilisant une representation du milieu par pixels, de 

nombreuses fonctions de transfert differentes appliquees a tout type de milieu heterosene 
peuvent etre facilement et rapidement culculees. 

La methode geometrique par exemple permet de trouver des dimensions de blocs 
equivalents permettant une tres bonne adequation du componement d'imbibition du bloc 
:o reel ou de la distribution des blocs reels, quelle que soit la forme (ou les formes) des 
blocs considered. La courbe de recuperation d'huile calculee sur la section de bloc 
equivalent, bien que simplifiee pur rapport aux methodes anterieures. se revele tres 
proche de celle calculee sur la section de bloc reel. 

D'autres caracteristiques et avantages de la methode selon ('invention apparaltront 
25 plus clairement a la lecture de la description ci-apres de modes de realisation donnes a 
titre d'exempics non limitatifs. en sc referant aux schemas figurant en annexe parmi 
ioqucls : 

- la figure I illustrc une' procedure connue permettant de determiner un milieu 
pre.sentant des fissures reyuliercs equivalent au milieu fissure reel. 
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- la figure 2 illu^tre la procedure >elon fin vent ion permettant de determiner un milieu 
presentant des fissures regulieres equivalent a un milieu fissure reel. 

- la figure 3 presence un exemple de pixels voisins servant au calcul de la \aleur 
attribute a un pixel. 

5 - la figure 4 presente un histogramme d'une distribution possible de pixels par rapport 

a la distance aux fissures. 

- la figure 5 presente une variation possible d'une zone envahie normee en fonetion de 
la distance aux fissures. 

- la figure 6 montre un autre exemple de pixels voisins dans trois plans differcnts Sk- 
10 1 . Sk et Sk+ 1 utilises pour un calcul 3D de valeurs attributes a un pixel. 

- la figure 7 presente une extension possible de pixels voisins permettant d'ameliorer le 
calcul des valeurs attributes a un pixel, et 

- la figure S montre. a des fins de validation, la bonne adequation entre deux courbes 
de recuperation d'huile OR(t) determinees en utilisant d'une part un bloc "en forme de 

15 peigne" reel et d'autre part un bloc rectangulaire equivalent. 

Une nouvelle procedure simplifiee permettant de calculer les dimensions de la section 
de bloc equivalente a la section "horizontale" d'un milieu fissure naturel est presentee ci- 
apres. 

En premier lieu, il convient d'indiquer que, selon 1'hypothese de fissures "verticales". 

20 e'est-a-dire perpendiculaires aux plans de stratification, le milieu matriciel est continu 
d'une couche geologique a une autre, et le probleme lie a la decouverte de dimensions de 
blocs equivalents devient un probleme bidimensionnel. Ainsi, le probleme auquel on est 
eontronte ici consiste a determiner la section carree ou rectangulaire equivalente de blocs 
matricieis numeriques pour chaque couche ou groupe de couches presemunt des 

25 proprietes de fissuration stmilaires. 

Deuxiernement. lequivalence d'un modele double porosite par rapport a un reservoir 
fissure doit etre etablie pour le comportement de I'ecoulemeni. Les ecouienienK au >em 
de reservoirs d'huile fissures sont esse ntie lie mem polyphasiques au cours de i'cxploitat ion 
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Ju iiisemem. avcc deux pnncipaux mecanismes de drainacc pour !a recuperation de 
I'huile dans la mat rice. ('imbibition captllaire et ie drainage par ^ravite. Les effets de ces 
deux mecanismes so conju-uent dans le cas de precedes de recuperation eau-huiic qui 
demeurent une strategic predominance dans le cadre du de'xeloppemcnt de nombreux 
reservoirs fissures. Des mecanismes compositionnels tels que la diffusion mtervienneni 
element dans les precedes de recuperation de gaz. Cest pourquoi la methodc 
geomeirique decrite ei-apres pour determiner les blocs equivalents est basee sur des 
concepts d ecoulemeni poK phasique. 

Un mode de realisation de i'in vent ion est decrit ci-dessous en relation avcc la rlcure 2. 
qui eonsiste sensiblement en la mise en adequation de la fonction de recuperation dhuiie 
R{t) (du milieu fissure nature!) obienue par la methode de reference citee et de la 
fonction de recuperation connue Req(t) pour le milieu transpose, pour un processus 
d'imbibition eau-huiie diphasique (au cours dun mecanisme de drainage par imbibition 
capillaire eau-huile). Cette mise en adequation est realisee pour chaque couche du milieu 
fissure, puis pour des ensembles de n couches. Dans ce cas. Ia fonction de recuperation 
R(t) obtenue est la somme des differences fonctions Rn(t) des n couches ponderees par 
les epaisseurs correspondanres Hn. Les fissures etant verticales. seules les dimensions 
horizontals du bloc equivalent sont determiners. La concordance des fonctions R(t) et 
Req(t) est done un probleme bidimensionnel. 

I Formulation geometrique 

Les fissures etant definies par les coordonnees de leurs points limites sur une coupe 
bidimensionneile XY dune couche. le processus d'imbibition par lequel de I'eau est 
presence dans les fissures et de I'huile est presente dans les blocs matriciels doit etre 
determine. On suppose que 1'invasion de la matrice par de I'eau est du type a piston. On 
admet que la fonction x=fVt) qui relie lavancee du front d'eau au temps est la meme pour 
tous les blocs matriciels. quelle que soit leur forme, et pour tous les blocs elementaires. 
Par consequent, la mise en adequation des fonctions Rft) et Reqit) equivaut a la mise en 
adequation des fonctions R(x) ei Reqfx). Ces fonctions defln.ssem physiquement des 
/ones normees envahie* par I'eau en fonction de lavancee du front d.mb^.t ion dans le 
milieu fissure. 
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En 2D. I 'expression anaKtique dc Reqjxj est la suivame : 
{Re^(.v) = I - 



a x h 



< i n 4 

(£/ - 2.v nh - 2x) - 2\ — + — \x -x~ , 



a b } n '■< h 



v ■= . 0. mini — - l 

L \2 2 J 



; Re<y( v ) = I. .v > mm 1 — . — ■ 



ou a et b sont les dimensions du bloc rectangulaire ou carre equivalent (a et b > 0) : 

La fo net ion Rt\) n'a pas dexpression analytique. EHe esi caleule'e a paiiir d'une 
discretisation de la coupe XY de la couche etudiee suivant I'algorithme defini ct-apres. 

2 Algorithme de calcul de la fonction R(x) 

La coupe XY de la couche etudiee est consideree comme une image dont chaque 
pixel represente un element de surface. Ces pixels sont espaces regulierement d'un pas dx 
dans ia direction X et dy dans la direction Y. L'algoruhme mis en oeuvre vise a 
determiner, pour chaque pixel de cette image, la distance minimale qui le separe de la 
fissure Ia plus proche. 

L'image est traduite par un tableau de nombres reels a deux dimensions : Pict[0 
:nx+I,0 :ny+! ] ou nx et ny sont les nombres de pixels de l'image dans les directions X et 
Y. En pratique, le nombre total de pixels (nx.ny) est par exemple de Tordre du million. 
Les valeurs des elements du tableau Pict sont les distances recherchees. 

Initialisation : Tous les pixels par lesquels passe une fissure sont a une distance nulle 
de la fissure la plus proche. Pour ces pixels, le tableau Pict est done initialise a la valeur 
0. Ceci est fait par un algorithme connu en soi (I'algorithme de Bresline par exemple) 
auquel on donne les coordonnees des pixels corrcspondams aux deux extremites d'une 
fissure consideree comme un segment de droite et qui initialise (h 0 dans le cas present) 
les pixels les plus proches. Les autres elements de Pict sont initialises a une valeur 
superieure a la plus grande distance existant entre deux pixels de I'imace. Cetie valeur est 
par exemple nx.dx+ny.dy. 

Calcul : Pour un pixel donne. le calcul de la distance recherchee a la fissure la plus 
proche se Ian u part ir des valeurs de distance deja calculees pour les pixels \oisinv On iui 
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alfccte une valour H ul m elle savere.jnferieure a la valeur qui lui a ete assignee 
initialenient. est le minimum des valeurs des pixels vo.sins auxquelles on ajoute la 
distance de ccs pixels a celui considere. 

Ce caicul esr realise en deux phases successives. Lors de la pas>e descendant, on 
parcoun l image ligne par Itgne. de haul en has et de gauche a droite ule Pict [ l.l ] p tc t 
[nx.ny] ). Les pixels dont on tient compte >om differents selon que Ton est dans une passe 
descendante ou une passe montante. Comme le montre la figure 3. les pixels en noir et en 
gnse.sont ceux que ion prend en compte respectivement durant les passes descendant 
et les passes ascendances, pour le pixel Px. 



Lecart oblique dxy etam defini comme : dxy = ^d\ 2 +dy 2 . Talgonthme s eerie 

pour j= I a ny 
pour i=l a n.x 

I Pict[i.j]= min Pictfi-I .j] + dx, : passe descendante 

i Pictfi- 1. j- 1] + dxy. 

i Pict[i.j-]J + dy, 

I Pict(i+ij-l] + dxy. 

I Pict[i,j] 
fin de boucle sur i 
fin de boucle surj 

pour j=ny a 1 , 
pour i=lx a 1 . 

I Pict[i.j] = minPict[i+lj] +dx. : passe descendante 
I Pict[i+i.j+l] + dxy, 

I Pict[i.j+1] +dy. 

I Pict[i-l.j+l] + dxy t 

I Pictfi.j] 
fin de boucle sur i 
fin de boucle sur j 

Histogramme : A pamr du tableau Pict ainsi calcule, on peut tracer un histogram™ 
en classant les valeurs non nulles (cclles affectees aux pixels hors des fissures; par ordre 
croissant, 

Le cumule de cet histogramme donne. pour toute distance delimitam deux .mervulles 
de rhistogramme. ie nombre de pixels non nuis dont la valeur est inferieure a cone 
distance. Dans Implication decnte a un milieu poreux fissure ou cette d.Mance 
correspond a I'avancce du Iront d'eau. le cumule de Thistogramme .ndique done lane 
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envahie pur 1'eau. La courhc Rtx) est obienue en divi.sani ce cumule par le nombre total 
de pixels non nuls (pour la normen. Lo nombre d'intervalles utilises cn abseisse pour 
I'histogramme correspond au nombre de points de discretisation de la courbe R{.\). On le 
choisit egal a 500 par exemple. 

3 Recherche des dimensions du bloc equivalent 

A ce stade, on connait la fonction R(x) et i'on cherche !es parametres (a.b) 
(dimensions du bloc equivalent qui minimisent la fonctionnelle i : 

N 

Jf a.b.) = £(R{Xj )-Req-a.b.xj )) 2 
i = l 

ou N est le nombre de points de discretisation de Rfx) et fx;) sont les abscisses de ces 
points de discretisation. 

Discretisation suivant I'ordonnee de R(x) : 

Pour donner autant de poids a tous les volumes d'huile recuperet lors de I'imbibition. 
la courbe R(x) est re-discretisee suivant un pas constant sur 1'axe des ordonnees (figure 
5). La suite (x;) utilisee par la fonctionnelle est deduite de cette discretisation. 

Minimisation de la fonctionnelle : 

Comme a et b jouent des roles symetriques dans I'expression Req(a.b.x). on utilise en 
fait la fonctionnelle : 

J(u.v) = Xl^f-O- Re?(u.v.jc ( )) : 



I R e q ( u. v.x) = u x x + v x .r : 

^ „ . so it i 

Re gin. v.x) < 1 



a b } 
_ -4 

a x b 



Mininuser cette fonctionnelle revient a trouver le couple (u,vj pour lequel 
J t u. V } = 0. Ceci est fait par un algorithmic de Newton. 
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Ensuitc. le coupio ta.b) recherche.est deduit de lu.vi. Tmis a.s peuxent so 
presenter : 

I ) v > 0 signifie quune de* valeurs du couple (a.b) est negative, ce qui na aucun 
>ens physique. On pose jlors v=0 dans I'expression de Req< li.v.x ). ce qui implique que 
les fissures sont parallels. L'operation est recommences et ie couple (a.b) est calcule 
com me suit : 



1 _ tt 
\h infm/ 



2) Le cas u* + 4 v < 0 est depourvu de tout sens physique egalement puisqu'il signifie 
que (a.b) ne sont pas reels. On pose alors u 2 + 4v = 0. ce qui impose que le bloc 
elementaire recherche a !a forme dun carre (a=b). Apres minimisation, le couple (a.b) 
e^t calcule comme suit : 

]a = « 

3) Pour les autres valeurs du couple ( u. v), on a : 



_ - « V m ~ + 4 v 
a — 



- - a - Vm : +■ 4v 
b = 



4 Validation de la determination d'un bloc equivalent 

La methode geometrique indiquee ci-avant pour une fonction de transfert d'imbibition 
a ete validee par rapport a une methode de reference convemionnelle d'un cout tres eleve 
qui necessite un simulateur decoulement polyphasique simple porosite discretisant les 
blocs matriciels et les fissures de telle maniere que les courbes de recuperation puis.seni 
etre comparees. Des simulations decoulement diphasique cnnvcmionnello uni c ;c 
realisees pour validcr les solutions fourn.es par la methode geometrique decrite plus 
haut. La validation peut comporter les eiapes suivantes : 



i: 
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to 



a i Caicul de la fonction de recuperation d'huile R(i) pour la section ueolouique (reellet 
avec la methode uon\ entionnelle I solution de reference): 

bt application de la methode geometrique a la section reelle. ce qui donne une solution 
\ a.b i: 

c) nou\elle utilisation de ia methode conventionnelle. caicul de la tone t ion de 
recuperation d'huile Reqm sur la section de bloc equivalent de dimensions (a.b) 
prealablement determinees. et comparaison avec la fonction de recuperation d'huile de 
reference R(w. 

La methode geometrique donne des dimensions de bloc equivalent permettant une 
bonne adequation avec le comportement a limbibition du bloc reel, quelle que soit Ia 
forme de bloc consideree. La courbe de recuperation d'huile calculee sur la section de 
bloc equivalent est toujours tres proche de celle calculee sur la section de bloc reel 
comme ie montre la figure 8. 

Une telle validation a egalement ete realisee avec succes pour une distribution de 
15 blocs presentant differences autres dimensions et formes. 

5 Autres applications de la methode 

Precision du caicul des distances aux fissures 

Dans l'algorithme de caicul des distances des pixels aux fissures a partir de I'imaae 2D. 
on peut ameliorer la precision du caicul en prenam en compte un plus grand nombre de 
:o voisins du pixel considered La figure 7 montre Ies pixels voisins done on tient compte 
quand on etend I'intluence d'un pixel Px a deux lignes et deux coionnes precedam ou 
sLiivant le pixel Px. De la meme facon, les pixels en noir et en grise sont ceux que I'on 
prend en compte respectivement durant les passes descendants et les passes 
ascendantes. ceux indiques par une croix etant elimines pour cause de redondance. 

:5 Pour augmenter encore !a precision des calculs. la zone d'influence des pixels pent 

Cure agrandie da vantage la 5 lignes et 3 coionnes ou plus). Cependant. pour I'utilisaUon 
presentee plus haut. une telle extension n'appone pas dans la pratique d'amelioration 
notable sur les resultats finaux. 
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Extension de la me t hod e a un objet tridimensionnel 

L'algorithme present e ei-dessus peut etre applique a un volume. Dans ce cas. chuque 
pixel represente un clement de \olume. Le tableau Pict est rempluce par un tableau a 3 
dimensions Pici/DfO :n\+I.O :n\ + l.O :nz+l] ou nx. n\ ei nz sont ies nombres de pi\eU 
5 en X. V ec Z. Pour le calcul en un pixel Px du plan horizontal nitmero k. Ies pixels vomiis 
dont on tient compte lors des passes descendante et ascendante sont represented sur la 
figure 6. 

Extension a une fonetion quelconque 

Dans I'exemple qui a ete developpe d'une etude dun phenomene dimbibition 
0 diphasique (eau-huile par exemple). on a cherche a determiner ia taille des blocs en 
relation avec les distances de points a la fissure la plus proche. La methode geometrique 
selon I'invention peut etre aussi utilisee pour d'autres types de transfers entre deux 
milieux contrasted comme par exemple les transfers thermiques entre un puits et un 
reservoir. Mais sunout. la fonetion "distance entre pixels" utilisee dans l'algorithme 
5 precedent peut etre remplacee par une fonetion quelconque reliant les points de I'image. 
II faut alors connattre, pour tout pixel de I'image. la valeur de cette fonetion entre ce 
pixel et ses voisins pris en compte pour le calcul. Cette fonetion peut. par exemple. 
exprimer les valeurs de transmissivites entre les mailles d'un reservoir dont les centres 
sont les pixels de I'image et servir au calcul du volume de drainage d'un puits fore dans ce 
0 reservoir. 

Dans un tel cas. Ies deux passes ascendante et descendante effectuees dans 
l'algorithme peuvent s'averer insuffisantes pour trouver une valeur minimale en tout pixel 
de I'image. On recommence alors jusqu'a ce que les valeurs calculees ne changent plus. 

En reprenant Ies notations presentees plus haut et en supposant que la fonetion 
5 F(i.j.k.l) renvoie la valeur de la fonetion entre Ies pixels (i.j> et (k.!>. l'algorithme 2D 
devient : 

change=vrai 
tant_quei change==vrai i 
change=faux 
o pour j= I a ny 

| pour i= I a nx 
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I | iemp=Pict[i.j] 

j | Picifi.j] = min (Pict[i-I.j) + Fii.j.i- 1 j). : passe descendant 

I I Ptct[i-I.j-1 | + F(i.j.i-|.j-| » 

! ! Pictfij-I] + F(t.j,i.j-1) 

i I Pict[i+l,j-l]-rF<i.j.i+l.j-l) 

I I Pictfi.il 

i I >i f Pictfi.jjotemp) change = vrai 

I tin boucle sur i 

t in boucle sur j 

pour j=ny a I . - 1 
t pour i=nx a I . - 1 
I j temp=Pict[i,j] 

! I Pict[i.j] = min (Pici[i+[.j] + Fti.j.i+I.j), : passe ascendante 

I I Pict[i+l.j+lJ + F(i.j.i+l.j+M 

I i Pict[i.j+I] + F(i.j.i.j+I ) 

I I Pict[i-I.j+IJ +F(i.j.i-1 j+l) 

! I Pictfi.j] - 

I [ si (Pict[i.j]otemp i change = vrai 

j tin boucle sur i 

tin boucle sur j 

fin cant_que 
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REVENDICATIONS 

I ) Methode permettant de simplifier la model isaiion d'un milieu geolouique poreux 
heterogene d'origine sous l;i forme dim milieu transpose ou equivalent al'in que le milieu 
transpose soil equivalent au milieu d'origine relati\ement a un type determine de 
tonet ion dc transfer! physique connu pour le milieu transpose, la methode comportani : 

- la formation dune image a au moins deux dimensions. (2D) du milieu geologique 
sous la forme d'une serie de pixels et I'association a chaque pixel dc la serie dune valeur 
mitiale particulate pour ladiie fonction. 

- la determination par etapes de la valeur a attribuer pour la fonction de iransfert 
physique a chaque pixel de ladite serie. en reference a des valeurs de la fonction 
attributes a des pixels voisins de 1'image. et 

- la determination dune propriete physique du milieu transpose ou equivalent par 
identification de valeurs de la fonction de transfer! connue pour le milieu transpose par 
rapport a la valeur de la fonction de transfert deierminee par etapes pour le milieu 
d'origine. 

2) Methode selon la revindication I. caracterisee en ce que le milieu geologique 
heterogene est traverse par un reseau irregulier de fissures toutes geometriquemcnt 
definies dans des blocs de formes et de tailles irregulieres. 

3) Methode selon 1'une des revindications I ou 2. caracterisee en ce que ladite 
fonction de transfert physique represente une distance entre differentes parties du milieu 
geologique. 

4) Methode selon 1'une des re vend icat ions 1 ou 2. caracterisee en ce que ladite 
fonction de transfert physique represente des transmissivitcs entre differentes parties du 
milieu geologique. 

5) Methode selon 1'une des revindications I ou 2. caracterisee en ce que ladite 
lonction de transfert physique represente des iranslens thermiques entre differentes 
parlies du milieu geologique lels que des traiisferts thermiques entre un reservoir el un 
puit.s uaversani ce reservoir. 
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o» Mcihode scion Tunc do revindications I ou :. caractensee cn co que laditc 
toncuon de transfer! physique represente tout transit de masse ou de flux entre 
dirferentes panies du milieu jieologtque. 

7) Methode selon la revindication 2, permertani de determiner, a partir d une ima-je 
dun milieu geoloizique poreux reel traverse par un reseau irreyuher de Insures, un milieu 
transpose comportant un ensemble de blocs reguiierement disposes et separe.s par un 
maillage regulier de fissures, ledit milieu transpose donnant sensiblement la meme 
recuperation (Req) de fluide lors dun processus d'imbibition capillaire que le milieu reel, 
la mcihode comportant : . • 

- la formation d'une image a au moins deux dimensions (2D) du milieu reel sous la 
forme d'une serie de pixels. 

- la determination, pour chaque pixel, de !a distance minimale separant le pixel de la 
fissure la plus proche. 

- la formation d'une distribution du nombre de pixels par rapport a la distance 
minimale au milieu fissure et la determination, a partir de cette distribution, de la fonction 
de recuperation (R) dudit ensemble de blocs, et 

- la determination de dimensions (a.b) des blocs reguliers equivalents du milieu 
transpose a partir de la recuperation (R) et de la recuperation (Req) du bloc equivalent 



2757957 





R(x)2Req(x) 
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a et b = dimensions correctes 
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